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The properties of the Hamiltonian developed in Paper II are studied showing that at a particular 
strain level a "localization" phase transition occurs characterized by the emergence of conjugate 
bands of coherently oriented cracks. The functional integration that yields the partition function 
is then performed analytically using an approximation that employs only a subset of states in the 
functional neighborhood surrounding the most probable states. Such integration establishes the free 
energy of the system, and upon taking the derivatives of the free energy, the localization transition 
is shown to be continuous and to be distinct from peak stress. When the bulk modulus of the grain 
material is large, localization always occurs in the softening regime following peak stress, while for 
sufficiently small bulk moduli and at sufficiently low confining pressure, the localization occurs in 
the hardening regime prior to peak stress. In the approach to localization, the stress-strain relation 
for the whole rock remains analytic, as is observed both in experimental data and in simpler models. 
The correlation function of the crack fields is also obtained. It has a correlation length characterizing 
the aspect ratio of the crack clusters that diverges as ^ ~ (sc — e)^^ at localization. 

PACS numbers: 46.50. +a, 46.65. +g, 62.20.Mk, 64.60.Pr 
Keywords: 



I. INTRODUCTION 

In Paper II of this series, we obtained the Hamil- 
tonian Ej{e,em) of a population of interacting cracks 
which is the energy necessary to lead a mesovolume of a 
disordered-solid system from uncracked and unstrained 
initial conditions, to a final crack state j at a maximum 
imposed strain Sm that is possibly different than the ac- 
tual strain e if the system has been subsequently un- 
loaded. Using this Hamiltonian, we prove here that at 
a well-defined strain Sc, the system undergoes a phase 
transition to bands of coherently oriented cracks. 

To study the nature of this localization transition, we 
must evaluate the partition function Z from which all 
physical properties depending on the crack distribution 
are obtained through differentiation. In Paper I, it was 
established that Z takes a standard form 

Z{e,e„.,T)=J2e-^^(^'^-y^ (1) 
j 

despite the fact that it derives from the initial quenched 
disorder in the grain-contact strengths and has nothing 
to do with fluctuations through time. The possible crack 
states j for a mesovolume are defined by a local order 
parameter ip{x) distributed at each cell a; of a regular 
square network of identical cells. The amplitude of (p{x) 
corresponds to the length of a local crack (always less 
than cell dimensions), and its sign indicates its orienta- 
tion (±45° relative to the principle-stress axis). 



Our approach for performing the sum over states be- 
gins by determining which fields (p maximize the Hamil- 
tonian. Because the temperature in strain-controlled ex- 
periments is negative, such maximizing states are the 
dominant terms in Eq. (1). Any change in the nature 
of the maximizing crack-fields or in the nature of the 
Hamiltonian in their neighborhood {e.g., the vanishing 
of a second derivative) corresponds to a phase transition. 

In Sec. II, the localization transition is identified and 
the geometrical nature of the crack fields in the "func- 
tional neighborhood" surrounding the maximizing states 
defined. In Sec. Ill, we sum only over this subset of all 
states to obtain an analytical approximation of Z. In Sec. 
IV, the free energy F = —TlnZ is differentiated with re- 
spect to e and T to determine both the sustained stress 
T, the energy U, and the entropy 5*. In the approach to 
localization, no singularities are present in either F or 
any of its derivatives with respect to strain or tempera- 
ture which demonstrates, among other things, that the 
stress/strain relation is analytic up to (and including) 
localization. In Sec. V, an external field J is introduced 
that couples to ip permitting an autocorrelation function 
to be obtained. All singularities at localization are in 
the second (and higher) derivatives of F with respect to 
J with the consequence that the correlation length di- 
verges as ^ ~ (sc — e)"^. 



II. PRINCIPLE OF THE TRANSITION 
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A. Extrema of the Hamiltonian 

We now determine the most probable states by max- 
imizing the Hamiltonian Ej{e,em) along the load path 



2 



e = £m ■ From the summary of Paper II, we have 

E, = E\e^) + (1 - g) [E^" {e„,m + S'"'(£,„)M} 

where is the energy of the intact material, E'^^ is 
the energy due to the crack field when crack interactions 
are neglected, and i?'"* is the energy due to crack in- 
teractions. The parameter q derives from the quenched- 
disorder distribution and is bounded as,l/2 <q <\. 

That the Hamiltonian must be maximized and not 
minimized comes from the temperature parameter being 
negative as was quantitatively established in Sec. IV of 
Paper I. Because we assume the system is intact before 
strain is applied, it is a fact of our model that the intact 
state is always the most probable. For this to hold true, 
the temperature must be negative in strain-controlled ex- 
periments because the arrival of cracks at constant strain 
always reduces the energy in a mesovolume. 

1. Mean-field terms 

A mean-field simplification of the model built in Paper 
II would reduce the Hamiltonian to the sole term 

+ (1 ~ q)E-- = i [aA2 + (1 - 

where A is the strain dilatation, 7 the shear strain, and 
a and combinations of the elastic moduli all as 

defined in Paper II. The second term is strictly negative 
and represents the weakening of the rock due to the crack 
porosity which is proportional to -0, the volume average 
of the positive field -0 — \ip\. Therefore this mean-field 
Hamiltonian is maximum when ip — which uniquely 
corresponds to the uniform intact state i/j = <y9 = 0. 

2. Interaction term 

The interesting term is the interaction energy E™^ . As 
defined in the summary of Paper II (the reader should 
consult this summary for the definition of all the terms 
in what follows), E™^ is a sum over wavenumbers k of 
orthogonal quadratic forms involving Rk and Ik, which 
are vectors containing the fc-space Fourier modes of the 
order-parameter fields (/? and ^. The sign of these forms 
is determined by the sign of the two eigenvalues of the 
symmetric matrices Pk- For any fc, at least one of the 
eigenvalues is positive, since [1,0] • Pk ■ [1,0]^ ~ Lk ^ 
A^K^ (1 - au|) > where 1/2 < a < 1 and Uk is a 
cosine. To determine the sign of the second eigenvalue, 
it is sufficient to take the determinant of Pk ■ Using -I- 

= 1, it is straightforward to show that 

dei\Pk\ = /^'^n\{l-a)[cvk+ojf ■ (2) 
This is strictly positive for every fc, except when 

Vk = sin(26lfc) = -uj/c (3) 



in which case the determinant and second eigenvalue 
are zero. The vanishing of the determinant is thus in- 
dependent of the norm of fc, and takes place at ei- 
ther of two conjugate angles 9^ — arcsin(— a;/c)/2 or 
9^ = 7r/2 — arcsin(— w/c)/2 where 9k represents the an- 
gle between k and the crack-orientation vector ei. The 
directions in fe-space at which the determinant vanishes 
will be denoted by the unit vectors fe*. Thus, the ma- 
trices Pk are positive definite; i.e., they have two strictly 
positive eigenvalues, except for those particular wavevec- 
tors lying along one of the two directions for which they 
become positive degenerate. The eigenvector of Pk asso- 
ciated with the zero eigenvalues is easily computed to be 
[1, -Mk/LkV. 

The positive-definite quadratic forms of i?'"* are multi- 
plied by a negative constant which implies that the max- 
imum of E™^ occurs when (^fc = 0fe = for every non- 
zero fe with the exception of those fe satisfying Eq. (|^). 
At these degenerate angles, the Fourier modes of (/? and 
-0 are related as 

Vk = -TT^fe- (4) 

Mk 

Now, the definition of the auxiliary field ijjx = \^Px\ im- 
poses a series of constraints between and i/'fc- The 
simplest is obtained by noting that the space-integrals of 
ip^ and -0^ must be the same which is equivalent to 

(jjki'-k - <pkv-kj = 0- (5) 

For a crack-state maximizing E™^ , this condition further 
requires that 

(^o-V^o)+ E i^-^^^k^-k = Q- (6) 
k = fcfc± " 

It will be seen momentarily that along the directions k"^, 
the factors l~-^^^±/^^^± ^-^e equal, and that this quan- 
tity is an increasing function of the shear-strain parame- 
ter Lu = (k37)/(Aki), starting at a strictly negative value 
when ui — (no shear deformation yet applied), and 
reaching at a particular value ujc- For every wavevector, 
^k'f-k = ||<^fe|| is trivially positive, and the definition of 
-0 also requires that (fg — ipQ < for any crack state. From 
Eq. (H), we can conclude that for uj < cUc, the only crack- 
states maximizing the interaction term E™^ must satisfy 
both (fil = 0Q and, for every non-zero k, ipk = '4'k = 0. 
Such a maximum thus corresponds to a spatially uniform 
crack field. 

At the degenerate point uj = cUc, the set of maximizing 
crack states goes through a drastic change. Any non- 
zero Fourier mode of ip and -0 along the directions k^ no 
longer modifies i?'"* so long as (^q — V'o! so long as 
the crack field has the same sign over the entire mesovol- 
ume. This degeneracy of i?'"* at w = Wc is at the origin 
of the localization phase transition. 



The critical value and the corresponding wavevec- 
tors fe for which non-zero Fourier modes of (/? and ip do 
not contribute to E™^ , are determined from the two con- 
ditions 



det(Ffc) = 
'l~Ml = 0. 



(7) 
(8) 



Using the solution of Eq. (|7|) given by Eq. in the 
definitions of and Mk given in the summary of Paper 
II, Eq. (0) then becomes an equation for uJc 



= 0. 



(9) 



From the definitions of Paper II, we have c > 1 while 
1/2 < a < 1. Thus, Eqs. can only be satisfied by 



= ±Vc2 - 1 



sin(20fc) = - Vc2-1 / 



(10) 
(11) 



With a radial confining pressure maintained constant, 
and a positive shear stress Taxiai > ^radial, the strain com- 
ponents of the rock satisfy £axial < ^radial and Eaxial < 
so that LO = (K3/Ki)(eaxial - eradial)/(eaxial + <^radial) IS a 

positive and monotonically increasing function of the ax- 
ial stress, until the rock possibly exhibits some positive 
volumetric strain (we will later show that this does not 
occur prior to localization) , where this quantity diverges 
to -l-oo and increases further startiiig from — oo. All of 
this establishes that Eqs. and ([g) have no solution 
until the first solution uj = a;+ is reached. At this par- 
ticular strain value, non-zero Fourier modes of Lp and "0 
having any wavevector lying in one of the two directions 
defined by Eq. (11) can be added to a mesovolume with 
no change in the interaction energy. 
For quartz as the rock mineral. 
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so that we find (eaxiai/CradiaOc — —1-2 at the transition. 
Our model thus predicts the localization transition to oc- 
cur after a sign reversal of ^radial but prior to the point 
where A = £axiai + ^radial changes sign. These results 
are consistent with what is observed in usual triaxial me- 
chanical experiments {e.g. ||^, ||, jsj). 

It can now be algebraically verified using the defini- 
tions of Lfe and Mf. given in Paper II, that l~-^f^±/^^.^± 
does not depend on the norm k nor on which of the two 
directions k"^ is selected. Further, it increases monoton- 
ically from a negative value to reach zero when ui = uj^ 
(facts used in obtaining the above results). 

B. Structure at the localization transition 

The goal here is to define the geometric nature of the 
states maximizing E™^ at the strain point uj^. Neces- 
sary conditions on the structure of the degenerate states 



axial 




FIG. 1: A part of the conjugate bands emerging at the crit- 
ical strain. The bands perpendicular to fc^ are exclusively 
composed of right-inclined cracks, while those perpendicular 
to fe~ contain only left-inclined cracks. 



were just given and these are easily made into sufficient 
conditions. First, the degenerate states must correspond 
to crack fields of constant sign. They thus satisfy ev- 
erywhere -0 — ifi or ip — ~ip or, equivalently, ipk — 'fik 
or V'fe = —'fk- Considering this together with the nec- 
essary conditions of Eqs. (^ and (0), requires that the 
degenerate states be one of two types: (1) (/? > ev- 
erywhere and the only possible non-zero Fourier modes 
of (fi have wavevector directions that satisfy Nk/Mk — 
Mk/Lk = — 1; or (2) (p < everywhere and the wavevec- 
tor directions satisfy Nk/Mk — Mk/Lk = +1. Using 
again the definitions of L^, M^, Nk given in the sum- 
mary of Paper II, the first type of degenerate mo de corr e- 
sponds to wavevectors satisfying sin(26'fc) = —\J(? — 1/c 
and cos(20fc) = —1/c, while the second type of mode 
has the same sine requirement, but an opposite value 
for the cosine. Using fe+ to represent the wavevector di- 
rection corresponding to the first condition, and k~ the 
wavevector direction for the second condition, we con- 
clude that the emergent degenerate crack states consist 
either of right-inclined cracks with spatial fiuctuations 
forming bands perpendicular to fc"*", or of left-inclined 
cracks forming bands perpendicular to k~ . Such geome- 
try is sketched in Fig. 0. 

These two sets of crack modes are conjugate to each 
other; i.e., symmetric to each other under inversion of the 
radial axis. Since they become statistically important as 
UJ ^ ujci whereas the intact state or uniform states are the 
important states prior to Wc, the system spontaneously 
breaks its symmetry at the transition, which is charac- 
teristic of a continuous phase transition. 

Further, the angle formed by these bands is at 45° — 
\(ik- \ from the axial direction. Using Eq. ( pi] ) and the 
definitions of ki and K2 in terms of the Lame parameters, 
it is found that this angle is typically between 15° and 
35° depending on the rock mineral ||j considered which 
is consistent with laboratory experiments. 
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Finally, we note that these special crack bands that 
leave i?"^* unchanged, make a negative contribution to 
the Hamiltonian through the mean-field energy E'^^ that 
is proportional to ip. Due to the range of elastic 
interactions, E™* is independent of the norm of fc (it de- 
pends only on its orientation). Thus, the spatial variation 
of the bands perpendicular to their lateral extent has no 
influence on i?'"* ; it only affects E'" through the number 
of cracks present. For large systems and a narrow band 
of only a few cell widths, tp — klj^ = kji where A and 
I are the linear size of a cell (grain) and of a mesovolume. 
Thus, such a thin band makes a negligible contribution 
to ip for large systems, and is energetically equivalent to 
the intact state. However, states with numerous and/or 
wide bands can make a non-negligible contribution to ip 
and are, therefore, less probable. So this transition in- 
deed corresponds to "localized" structures. Only those 
states with a small number of small width bands along 
the special directions are the statistically emergent ones 
as is observed in actual experiments on rocks. 



III. OBTAINING THE PARTITION FUNCTION 

The sum over crack states in Eq. (|l|) is equivalent to 
the functional integration 



-i5[y,£,£„]/T 



(12) 



Since our Hamiltonian is expressed in terms of the Fourier 
modes (pk, it is shown in standard textbooks |H, |^ that 
Z further transforms to 



(13) 



fcex 



where and are the real and imaginary part of ^fc, 
and T is a half-space of the set of the wavevectors cor- 
responding to the non-zero modes; i.e., corresponding in 
2D to the discrete set (fci,/c2) — (^■nlln\,2'Elln2) with 
(ni,n2) e 'S? . There is a small- wavelength cutoff given 
by max(|ni|; |n2|) < l/K that ensures that Lp does not 
vary on scales smaller than that of a cell, and there is 
the arbitrary criterion fci > made to divide this space 
into two symmetrical parts. Equation ( [T^ ) is valid up 
to a multiplicative constant that has no physical impor- 
tance since the properties of a system correspond to the 
derivatives of the free energy F = —TlnZ. 

An analytic approximation for Z is obtained by per- 
forming the functional integration over a properly chosen 
subset of all the possible crack-states. The definition of 
this subset is based on what was learned in the preceed- 
ing section; namely, that among the states having a given 
non-zero crack occupation ip, the most probable are the 
uniform states, and precisely at the phase transition, cer- 
tain banded states may arrive at almost no energy cost, 
and these emergent states also have the same sign over 



space. Thus, the geometrical characteristic of all such 
states in the "functional neighborhood" of the minimiz- 
ing state is that in each one, all cracks are oriented in the 
same direction (either left or right). This property justi- 
fies making a so-called "constant-sign" (or "mean-phase") 
approximation for the partition function in which only 
those states in which the sign does not change in space 
will be considered. This still includes a huge range of 
states in which (p spatially varies. The excluded states in 
this approximation are guaranteed to have lower proba- 
bilities than the included ones and, as such, should have 
a negligible influence on the physical properties of the 
system. In this aproximation, the Fourier modes of the 
auxiliary ■0 field are trivially related to those of f as ei- 
ther ipk — 'Pk for the positive states, or V'fc — —'Pk for 
the negative states. 

We now rescale the temperature as T = h^T' /t^ . 
From the definition T = dU/dS and the fact that U is 
an energy density independent of I while S is extensive 
and thus increases as , we have that T scales as . 
In taking the thermodynamic limit in what follows, it is 
convenient to work with the purely intensive parameter 
T' (that is independent of £ ). Our partition function 
within the constant-sign approximation then takes the 
form 







Zc:ijv(pexp< 








+ j Vip exp < 





d+e 



d+e 



'Po_ 

I'D 



fceT 



fcGT 



p>k 

I'D 



>(14) 



where V^p is a compact notation for the functional mea- 
sure d(pQ rifcexl'^'^fe'^'^fe)! ^iid where and /_ represent 
integration over the subsets of (/? fields that are every- 
where either positive or negative. The quantities d, e 
and are defined in the summary of Paper II as 



aA^ + (1 - a)7^ 



(15) 



/I 

\{Lk±1Mk + Nk){e) 



-q{Lk±2Mk + Nk) {£„,)]. (17) 



Recall that the values of the actual strain e intervening 
in the probability distribution and in the partition func- 
tion are those along the load curve for which e = 
Their formal distinction only plays a role when par- 
tial derivatives of the free energy are taken to define 
stress. We note then that the value of at e = Em 
is w± = -(l-g)e2[l,±l]-P(fc)-[l,±l]^/(l-a) and since 
we have shown that F is a positive-definite matrix, and 
that the temperature T' is negative, we have that w'^ jT' 
in Eq. ( p^ is strictly positive. 

The symmetry of the problem under the parity trans- 
formation (inversion of the radial axis) guarantees that 
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both integrals in Eq. ( p^ are equal. Accordingly, only 
the first integral over positive crack states will be treated. 
This integral separates into products of Gaussian in- 
tegrals with the only remaining coupling between the 
Fourier modes coming from the complicated constraints 
on the integration domain boundaries that are what guar- 
antee to have the same sign everywhere in real space, 
and -0 to lie within [0,1]. But in order to study any 
singular behavior of the free energy F in the vicinity of 
localization, Z is determined in the thermodynamic limit 
in which both the system size and mesovolume size i are 
taken to be infinite. In this limit, the complicated inte- 
gration bounds in fc-space are not relevant. The integra- 
tion can be carried out entirely on K+ for i^q/I^ , and M 
for each of the variables ip§/£^, Vk/^^ without chang- 
ing the result because the contribution to these integrals 
in the thermodynamic limit comes from the immediate 
neighborhood of ^o/^'^ = and ^kl^^ — 0. 

A technical proof of this can be obtained as follows: 
using R"^ and R as the integration domains produces 
an upper bound for Z since this includes every posi- 
tive crack-field. A lower bound can be obtained by re- 
ducing the integration domain to a subset of the set 
of all positive crack-fields in which < (^q < and 
Efcex (l^fe I + l^il) < min[(^o/V2, i^P - ^o)/V2]. Inte- 
grating mode by mode over this polyhedra, the result can 
be shown to be asymptotically equivalent to the result of 
the upper bound in the limit where becomes infinite. 
This exercise is left to the attention of the reader. 

Thus, no coupling between the fc modes exists in the 
thermodynamic limit, and our approximation of the par- 
tition function takes the convenient form 



Z^2zoZiY[ [z(fc)2] 



feST 



where zq = e ^ '^^^ ^ and 



21 = 



2;eR+ 



(18) 



(19) 



(20) 

In the Hmit I -l-oo, these two integrals become 

zi - h°T'l{t^e) (21) 

z(fc) - ^T^korj^' l\i^w+{k)\. (22) 

Using Eq. (|l^), one then obtains the free-energy density 
in the thermodynamic limit 

F = -T'{\nZ)K^ /(.^ 



ID 



ker 



In 



A^T' 



(23) 



The contribution zi has vanished in this limit due to 
the fact that a;lnx^Oasa;^0. This is a technical 
consequence of the fact that for states composed of a few 
single bands, ip vanishes in the thermodynamic limit, as 
commented upon in the previous section. 



IV. SYSTEM PROPERTIES AT LOCALIZATION 

The remaining task is to link this free energy to the ob- 
servables of the system by taking the partial derivatives 
of F in the limit as localization is approached. 

The two partial derivatives of primary interest are 
those that give the dimensionless entropy density s — 
A^S/iD and the stress t. From Paper I, we have 



dF 



and T = -— 



dF 

de 



T'e„ 



The free energy of Eq. ( p3| ) is rewritten by replacing the 
sum over the wavevectors X^fcex with a continuous in- 



tegral ^^/(27r)-° kdkJ^dO. After performing the 
trivial integration over dk we have 



where / is the integral 



/ — TT In 



A^r 



1= I \n{-w+)de. 

-■it/2 



(24) 



(25) 



The integrand is a temperature-independent strain 
function so that —dF/dT' gives 



/ TT 

^='2+2 
while from F — U — T's 



1 + ln 



U = d+^T'. 



(26) 



(27) 



Since d represents the linear elastic response of an in- 
tact rock, and T' decreases from zero to negative values 
as damage accumulates, this expression shows that the 
average energy decreases due to the presence of cracks 
and is thus consistent with the negative curvature of the 
strain/stress load curve observed experimentally. 

Before addressing how s and F (and their derivatives) 
behave at locaHzation, we first establish the stress and 
temperature behavior at localization. 



A. Mechanical behavior at localization 

Consider the stress components a — —2dF/d^ and 
p = —2dF/dA where a (shear stress) and p (pressure) 
are both positive and related to the axial and radial stress 
components as 



and 



P 



(28) 



In standard lab experiments, the axial stress varies 
while the radial stress = —pr is kept constant. The 
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strain components 7 (shear strain) and A (dilatation) are 
similarly related to the axial and radial strain as 



7 



and Is. — Sa 



(29) 



Using the definition of w+ [Eq. (^^] along with the def- 
initions of Lfc, Mfc and Nk given in the summary of 
Paper II, we differentiate the integral / with respect to 
the actual strain variables, evaluate along the load path 
(Am = A and 7™ = 7), use the definition uj — [J37/A 
with the new constant 073 — ks/ki and make the change 
of integration variables z = tan"^ 6 to obtain exactly 



L03 



{l-q)A 
1 

(I-'Z)A 



{2ti - ujd^I) 



(30) 
(31) 



where q = 1 — l/(fc + 2) is the constant associated with 
the exponent fc > of the quenched disorder distribution, 
and the integral d^^I is defined 



duiQ dz 



(32) 



g 1 + z2 

with g{bj, z) given by 

g{uJ,z) = [l-a-2(l-a)c+(l-a)c2+w2]z'* 
+ [Aauj + 4(1 - a)cuj]z^ 
+ [2 + 2a + 2(1 - a)c^ - 2{2a - l)uj^]z^ 
+ [-4q;lj + 4(1 - a)cu]z 
+1 - a + 2(1 - a)c + (1 - a)c^ + lo^ . (33) 

Thus, the shear stress and pressure can be written 

r 



-2(1 - a)7 



CJ3 



A(l-q) 



dul{u) (34) 



2aA 



r 



(l-g)A 



[27r - Lod^Iiuj)] . (35) 



The integral duil is solved using the residue theorem once 
the roots z of the quartic g{uj, z) have been found. 
This quartic decomposes into the exact form 

g{uj,z) [z - C(w)] [z - C {^)]u{uj, z) (36) 
u{u:,z) - p{io)[z~i{Lo)][z~C{oj)] (37) 

where the star indicates taking the complex conjugate. 
The roots ^(1^^) and C*('^) both merge to the real axis in 
the approach to localization w ^ Wc, while the other two 
roots and £,*{uj) remain complex at localization. 

There are thus three simple poles C,{uj), ^(w) and i 
contributing to duji if the loop is closed in the upper-half 
z plane so that the residue theorem yields 



5.5(C) 



+ 



(38) 



where Q{ } designates taking the imaginary part. We are 
interested in evaluating this integral (and therefore the 
roots C and ^ and the function p) only in the approach 
to localization; i.e., when 5uj = lo — ujc can be considered 
small. In this limit, the second and third terms of Eq. 
(^sl) (the residues from ^ and i) have numerators and de- 
nominators that are both order in 5u so that it suffices 
to know the behavior 



(39) 
(40) 



However, the residue related to C is proportional to Slu 
in both the numerator and denominator which requires 
knowledge of this root to second order 



C(w) = Co + Ci'Jw + C25lo^ 



(41) 



The various strain-independent constants ^i, pi, and Ci 
are all known groupings of the elastic constants derived 
from Eqs. (|3|), (||) and (|3^). The final result for the 
integral after an enormous algebraic reduction is 

d^I = Ic + hSoj (42) 

where the constants Ic and Ii are exactly 

2-c2 



2tt 



and Ii = 2tt- 



(43) 



1. Stress and strain at localization 

The shear stress and pressure may be written 
cr = (To + cr'"' and p ^ po + p'"* 

where ctq = —2(1 — a)j and po = —2aA are the trivial 
linear variations of the uncracked material. We have just 
shown that at localization {Suj = 0), the non-trivial shear 
stress due to cracks and crack interaction is 



_int 

CT^ = — 



(1 - g)A, c2 



< 



while the non-trivial pressure is 
:„t 2ttT' 



Pc 



(1 - g)A,c2 



< 0. 



(44) 



(45) 



That these critical values are both negative follows be- 
cause (scaled temperature at locaHzation) is nega- 
tive and Ac (total dilatation at locaHzation) will soon 
be shown to be negative. Equations (^J) and ( |45| ) say 
that the presence of cracks has lowered both shear stress 
and pressure relative to an intact material at the same 
strain. This is indeed what is observed in experiments. 

To quantify the nature of Ac, we use that the confin- 
ing pressure pr is a known positive constant in standard 
experiments on rocks so that 



Pr 



-aAc + (1 - a)"fc 



T^[2Tr-{uj, + ujs)I, 
2(l-g)Ac 



(46) 
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Together with ujc = bj^^d ^c, this represents an equation 
for Ac 



^3 



T'^ [27r - [u, + cj3)/c 
(1-9) 



(47) 

Because T' varies with strain, we have that T'^ is also 
a function of Ac so that Eq. ( ^7|) is more than a sim- 
ple quadratic in Ac. To obtain an order-of-magnitude 
estimate of T'^, we use the approximate temperature ex- 
pression based on non-interacting cracks 



1 

tI 



X In • 



2A2 



d2„(l-g)[K2+K3(WcM)2]A? 

2r 



(A + 2^)dmA2[K2 + K3(wc/W3)^ 



19/(1-9) ^ 

-1 (48) 



After putting Eq. (48) into Eq. (p7|), Ac is numerically 
determined using Newton's method. The predicted Ac is 
negative for the range of confining pressure Pr of inter- 
est and remains negative for all ranges of elastic moduli 
found in rocks. The signs of the various terms in Eq. 
( ^ imply that the transition happens when the tem- 
perature has sufficiently departed from zero, but is still 
negative. Typical results from the numerical evaluation 
are ~ — 10~^(A -I- 2/x), which confirms the rough esti- 
mate given in Sec. V of Paper II. The typical value for 
Ac is a few percent; i.e., the order of magnitude experi- 
mentally observed at peak stress. 



The conclusion is that at localization, both dilatation 
Ac and shear strain 7c = tOc^d^z are negative while 
1 7c I » |Ac|. This demonstrates that the radial strain 
= Ac — 7c is positive at localization, which is also 
consistent with experimental observations. 



2. Stress, strain, and temperature derivatives at localization 

We now address how the stress and strain components 
as well as the temperature are changing with the negative 
of axial strain e ~ —Ea = — (A + 7)/2 at localization. 

In the approach to locaHzation we write A = Ac -I- 5A, 
7 = 7c (S7, and T' = T'^ + ST' using the exact dif- 
ferential equation for temperature to define ST' in what 
follows (not the approximation). The condition that Pr 
is constant requires that 



SAl-a 



T^ [2tT - {2lUc + LU3)Ic - (Wc + 



-f57|l-a + 
ST' 



2(l-g)A2 

rc^3 [/c -I- (uJc + CJ3)/i] 

2(l-g)A2 

{277 - (Wc ^ W3)4} = 



2(l-(7)Ac 
which along with —ISe = 5A + S^ gives 



IcfA ^ 2(1 - g)(l - a)Ag/rc^ + ujj, \h + (c^c + + ^di^T'd [2ti - {uj^ + c^3)/c] dT'/de 

2 de -2(1 - g) A2/T^ + 2^ - 2{u, + uj3)Ic - {uJc + ^3)2/1 

1^7 ^ 2(1 - q)aAl/T^ -271+ {2oj, + uj3)Ic + (^c + ^3)^c/i - Ac/(2r^0 [2^ - {uj, + U3)Ic\ dT'/de 

2de ^ -2{l~q)AdT^ + 2TT-2{LUc + uj3)Ic-{LOc + uj3yii 



(49) 
(50) 



To obtain an exact expression for dT'/de (within the context of having employed the mean-phase approximation), we 
use the formaHsm of Sec. IV. A of Paper I to write 



dT'U+ (OaU + dA,„U + I) 9t'A + (a^J7 + d^^U+ |) dr- 



dT' 

de 



(51) 



Using Eq. 



for U, we have Ot'U = 7r/2, OaU = aA = -po/2, d^U = (l-a)7 = -(Jo/2, Oa^U = 0, and d^^U = 



so that the temperature derivative at localization is given by 



[(27r - u;Jc){l -a)+ LUsI^a] 2(1 - q)Ac + usiu;^ + t^3)(27r/i + LOj.PdT'J A 



IdT' _ 

2 ^ " -[27r - (cjc + cJ3)/c]2 + (1 - g)7r [-2(1 - q)Al/T'^ + 2tt - 2{uj, + cj3)/c - {^c + ^zfh\ ' 

I 



(52) 



This derivative is numerically calculated to be finite and 
negative for the ranges of elastic moduli and radial con- 
fining pressures of interest, thus indicating that the lo- 
calization transition always preceeds the phase transi- 



tion where the temperature diverges to —00. Since rocks 
fail immediately after localization, the temperature- 
divergence transition is not observed in rock experiments. 

Last, we determine the variation of the stress compo- 
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B. Entropy and its derivatives at localization 

The exact result d^^I = Ic + hSuj with Ic and /i as 
given by Eq. means that the integral / of Eq. ( p5| ) is 
itself both finite and continuous in the Hmit as Suj ^ 0. 
Because it has further been shown that T' remains finite 
and continuous at localization, Eqs. ( |24| ) and (p6|) then 
show that both the free energy and the entropy (and all of 
their derivatives with respect to strain) remain finite and 
continuous as Slj — > 0. This demonstrates exactly that 
the locaHzation transition is a continuous phase transi- 
tion and allows us to classify it as a critical point. 



CORRELATION FUNCTION 



A. Derivation of a diverging correlation length 



FIG. 2: The localization value of the axial tangent modulus 
da/de as a function of the radial confining pressure Pr. The 
three curves represent different assumed bulk moduli for the 
mineral. The other rock properties are F = 10 J/m^, dm ~ 
10/im, /I = 15 GPa, and q = 3/4. 



nents with axial strain e at localization. Since Pr is con- 
stant, we have that dp/de — da/de = —dra/de. These 
derivatives define the so-called "tangent modulus" given 
by 



1 da 

2de 



W3 



2(1 - q) Ac de 



2(1 



(1-a) 



T' 



2(1-<7)A 



2 '^3 



g)A2 

d'y 
de 



de 



^ (53) 



where the derivatives dA/de, dj/de, and dT' /de have 
been given above. 

In Fig. H, we plot how da/de varies with radial con- 
fining pressure for various values of the elastic constants. 
The plot shows that for a sufficiently large ratio of bulk 
to shear modulus, the axial pressure is always decreasing 
at localization, which means that it has already passed 
through the stress maximum. However, for sufficiently 
small bulk moduli and at low confining pressures, lo- 
calization can also occur prior to peak stress. Thus, 
peak stress and localization are distinct in our theory. 
LocaHzation can occur in either the hardening or soft- 
ening regime depending on the bulk modulus and con- 
fining pressure. When localization occurs in the soften- 
ing regime (large bulk modulus), the strain/stress curve 
around peak stress is necessarily an analytic (quadratic) 
function, whereas when it occurs in the hardening regime 
(small bulk moldulus with small confining pressure), the 
peak stress presumably corresponds to a sharper varia- 
tion as microcracks start to coalesce along a weakened 
band and unstable failure sets in. These predictions are 
consistent with the experimental observations. 



The qualitative study of Sec. pB| lead to the conclu- 
sion that the localization transition is associated with 
the creation of conjugate bands of coherently oriented 
cracks. In this final section, the statistical correlation 
between cracks will be quantitatively addressed. 

The autocorrelation function is defined 

G{x, v) - {ip{x)^{v)) - {^{x)) My)) (54) 

and will be determined using a standard method of sta- 
tistical mechanics [|[ js). First, the Hamiltonian E[ip] 
is generalized to include an additional coupling of the lo- 
cal field ip{x) with an aribitrary field J{x) coming from 
some external source 

E'[ip,J]=E[ip]- I d^xJ{x)Lp{x). (55) 

The partition function becomes then a functional of the 
external field 



z[j] = / n (^^^ 



(56) 



and the averages involved in Eq. (|5j) are obtained by 
taking functional derivatives of Z[J] with respect to J 
and then letting the external field go to zero; i.e., 



{(p{x)) = lim 



T SZ 



{ip{x)ip{y)) = lim 



7^0 Z SJ{x) 

S'Z 



/^o Z 6J{x)SJ{y) ' 



(57) 
(58) 



Since the original Hamiltonian is most easily handled in 
Fourier form, the external coupHng will be expressed as 



d^x J{x)if{x) 



-^(jo^o + 2j2 Jk^k + 2 Jk^i ] (59) 
\ ker fceT / 
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where the superscripts R and / refer once again to the 
real and imaginary part of a complex quantity. The func- 
tional derivatives relative to J{x) must then be expressed 
by their counterparts in Fourier space 



S 



5Jk 9 



SJi d 



5J{x) \5J{x) dJl^ ' 6J{x)dJ{ 

^ ' fceTu{o}\ ^ ' '-''^k ^ ''-''^k, 



^ ( . d .... a 



feeTu{o}\ 



cos(fc-a;) 



sin(fc-a;)-^ 1.(60) 



The modified partition function will again be determined 
using the constant-sign approximation, but now the pres- 
ence of the external field breaks the symmetry between 
the sum over positive and negative crack-fields, so that 
both terms need to be kept in the generalization of Eq. 
(|lj). This leads to a slightly more complicated version 
of Eq. (jl|) for the expression of Z in the thermodynamic 
limit 



Z 



[z+{k)zt{k)\+z^\{ [z^(fc)z7(fc)] 

I fceT fc6T 



where zq is again the trivial intact term, and where 



(61) 



z%{k) = 



a:eR+ 

dx e 



2J«a:+£"u;±(fc)a:21/A"T' 



(63) 



with zf{k) having the same form as z^{k) after replac- 
ing with J^. In the following, the forms implying 
derivatives with respect to are to be implicitly un- 
derstood as having the same form as their counterpart 
with respect to (this imaginary components will not 
be explicitly written out). 

The integrals are easily performed giving 



A^T'/ 



ztik) 



exp 



e^^ ± Jo 



£Dw^{k) 



(64) 



(65) 



The first derivatives of Z with respect to the external 
field are then 



I — ^""^1 z+ n [z+{k)z+{k)] 



dJo 



A^r 



{e£D - Jo)' : 



lU [^«(fc)z7(fc)] (66) 



ex 



= za< z{ 



ker 



2 1 

AW(fc)T- [^^W^7(fe)]|(fi7) 



Letting the external field go to zero, both of these terms 
disappear, so that using the chain rule of Eq. (|60|), the 
average of the crack variable >p at any point a; in a meso- 
volume is given by Eq. (|5^ ) to be 

{^{x)) = 0. 

As expected, there is no spontaneous symmetry breaking 
prior to the transition. 

Consequently, the autocorrelation function reduces to 
only the second derivatives of Z in Eq. (^). Differenti- 
ating Eqs. (|66|)-(67) with respect to Jq, J^i and J^, , and 
taking the limit where J goes uniformly to zero leads to 



92 Z 



dJkdJk' 



d'Z 



d'z 

dJidJl 
2A°T' „ 



1 



1 



-(fc) 



(fc) 



kk' 



AD£DT' 



^0 



where Zq — Z[0] is the original partition function without 
external source. All the remaining cross-derivatives go to 
zero 



d'Z 



d'^Z 



d'^z 



dJ^dJo dJ^dJo 



0. 



Through the chain rules of Eq. ( pOD , these equalities show 
that the autocorrelation function has the form G{x; y) = 
G{x — y) due to the symmetry of the problem under 
translation for an infinite system. The Fourier transform 
G{x; y) = Y.k Gfce*'^ (^-3'V£^ is thus given by 



Gfc = 2A'°r' 



1 



1 



-(fc) 



ik) 



(68) 



when fe ^ 0. The special value Go = 2A^^T'^ /e^£^° 
does not play any role in the thermodynamic limit. 

In real space, the autocorrelation function is obtained 
by an inverse Fourier transform 



Gix) = 



1 

4^ 



27r/A 



kdk 



d9G{k)e 



ikx 



(69) 



Using G{—k) = G{k) which is a consequence of w''' being 
TT-periodic functions, and working in polar-coordinates 
X — [x, Ox) and fc — (fc, 9), the angular integral is divided 
into two symmetric domains which gives 



'^/^kdk 
2-k'- 



G{x,ex)^l —I d9G{k,e)cos[kxcos{e-9x)] 



Since and therefore G only depend on the angular 
part of fe, the integral over fc = |fe| yields 



Gix,9x) 



-^^^G(.) 



27r2 



2n 



■ sm 



/ 2ttx cos{9 — 9x) 
Aa;cos(6' - O.^) \ A 

1 f / 2Trx cos{9 — 9x) 







COS 



A 



- 1 
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For X ^ A, this integral is dominated by a neighborhood 
of9 = 9x + n/2, of angular size ciA/x with ci a constant 
of order unity. The function G{d) is almost constant 
over this small neighborhood, and this integral can be 
well approximated as 

G{x,e,) = — Ig [j) 

with the dimensionless integral Ig defined as 

27r sin (27ru cos{9)) 



Ig{u) = 27r^ / dO 

{cos (27rwcos(0)) — 1} 



An asymptotic study of this oscillating integral for u ^ 1 
shows that Ig{u) ~ C2/M, with C2 a positive constant of 
order unity. Reformulating Eq. (pSh with 



h{e) = 2C2T' 



1 



1 



w+{e + n/2) w-{e + T:/2) 



(70) 



gives the real space autocorrelation function in the form 



G{x,0.,) ^ h{e,) 



A 




(71) 



FIG. 3: Form of an iso-autocorrelation curve in the approach 
to the localization transition. 



the branches remains finite; i.e., that the aspect ratio of 
the branches £,/p also diverges as {ec — e)~^. This is qual- 
itatively illustrated in Figure 3. This prediction can be 
interpreted as corresponding to the formation of clusters 
of microcracks having aspect ratios £,/p that diverge as 
the cracks organize into long thin structures along which 
the sample will ultimately fail to form the experimentally 
observed shear bands. 



This establishes that along any direction, the autocor- 
relations decay as A/x (for two points separated by a 
significant number of grains, x :s> A). 

Concerning the angular dependence of G, the sym- 
metry of the system under parity leads to w~{9) — 
w+(7r/2 — 6*) (which can also be verified directly from the 
definitions of and the dependancies of Lk, M^, Nk on 
Uk = cos(26') and Vk = sin(26') given in paper II, together 
with the fact that the parity symmetry keeps v constant 
but changes the sign of u). This, along with the tt peri- 
odicity of w"^, shows that G is symmetric under parity; 
i.e., G{x,7T/2-d^)^G{x,e^). 

The angular dependence is best shown by considering 
curves of iso-correlations G{x,9x) = C3 where C3 is con- 
stant along a curve. Such curves obey x — Ah{9x)/c^. 
The direct study of the function w'^ shows that it admits 
quadratic maxima along the directions ^'*'[7r], scaling as 
max(w+) — w+(0+[7r]) = -~a{5Lo)'^ when the transition is 
approached, where a is a positive constant. This comes 
from the fact that E™^ is degenerate exclusively for the 
critical angles 9^, at reduced strain ujc. Outside a small 
neighborhood of 0"'"[7r], «;+ remains bounded. The defi- 
nition of h and the exchange under parity of w"^ shows 
then that such an iso-correlation curve has four branches 
(spikes) along the directions ±0+ ± 7r/2, whose extent ^ 
diverge to +00 as 

C ~ 2A(-r') — (c^c - u:)'^ ~ C4A(£, - e)-\ (72) 

The fact that remains bounded outside any small 
neighborhood of 9'^['k\ also means that the width p of 



B. Experimental measurement of ^ 

It would certainly be desirable to have direct experi- 
mental verification of whether the crack bands have as- 
pect ratio that diverge as l/(£c — e)^- Unfortunately, 
there are many practical problems that have prevented 
the direct measurement of the autocorrelation function 
of cracks in materials like rocks. We comment here on 
three types of measurements that either have or could be 
used to quantify the autocorrelation. 

First, following ideas used by Davy and Bonnet |Q in 
interpreting their sandbox shear experiments, one can 
measure the local deformation of a large sample by cov- 
ering the surface with pixels and mentoring the shear 
strain of each pixel. The total shear strain of the sys- 
tem is then approximated by taking the average over the 
surface pixels. If the system deformation is plotted as 
a function of the pixel size, it is expected that when the 
pixels are smaller than the emergent band structures, the 
system deformation will decrease as a powerlaw of in- 
creasing pixel size as was observed by Bonnet and Davy. 
However, at a particular pixel size there is a cross over 
to a constant system deformation as pixel size increases. 
The pixel size at the cross-over point is at least an indi- 
rect measurement of the correlation length ^ above which 
a volume-averaged description of the system holds with 
properties independent of the pixel size. 

Second, a direct measurement of the autocorrelation 
between cracks can in principal be obtained via acoustic- 
emissions monitoring |10|. However, the present resolu- 
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tion of this method (minimotors in contimotor-scalc spec- 
imens) and the difficulty in determing the mode of the in- 
dividual crack events prevents having a satisfactory sam- 
pling for statistical analysis. It seems that improvements 
on these present limitations are possible. 

Last, by analogy with the probing of spin pop- 
ulations by electromagnetic waves to study the 
ferro/paramagnetic transitions, it should be possible to 
send plane sound waves through a system and measure 
the scattering cross-section as the waves scatter from the 
structure of the evolving microcrack population. We have 
not yet obtained the rigorous connection between such a 
measured cross-section and the Fourier-transform of our 
autocorrelation function; however, such a relation almost 
certainly exists. No experimental attempts to measure 
the correlation function of cracking systems in this man- 
ner has been attempted to our knowledge. 

VI. CONCLUSION 

We now summarize the principal results that have 
emerged in our study. First, we have demonstrated that 
at a well-defined strain point w = Wc, thin bands of co- 
herently oriented cracks can be added to the system at 
no energetic cost. Such locahzed structures break the 
symmetry that held when w < Wc and correspond to a 
phase transition that we named the "localization transi- 
tion". It was demonstrated that the free energy F and 
entropy of the system remain continuous and finite at the 
localization transition which justifies caUing it a critical- 
point phenomena. Such continuity also demonstrates 
that the stress/strain behaviour of the rock is entirely 
analytic up to and including localization. The only di- 
vergence at localization is in the second derivatives of F 
with respect to the external field J. The consequence is 
that the correlation length (aspect ratio) of the emergent 
crack clusters diverges as (wc — Presumably, if the 



"mean-phase" approximation had not been invoked and 
if order-parameter contributions proportional to ip'^ and 
higher had been retained in the Hamiltonian through a 
renormalization scheme, then a non-trivial exponent on 
this scaling law might emerge. 

The mechanical behavior of the system at locahzation 
exhibits many qualities observed in actual experiments 
on rocks. First, the stress components at localization 
are reduced relative to their values if the rock had re- 
mained intact. The total dilatation Ac remains nega- 
tive at locahzation, even though the radial strain is posi- 
tive. With radial confining stress kept constant, the tan- 
gent moduli da/de are, most normally, negative at local- 
ization indicating that the load curve has already gone 
through a smooth quadratic peak stress prior to locahza- 
tion. Nonetheless, for rocks with a sufficiently low bulk 
modulus and at sufficiently low confining pressures, the 
localization can occur in the hardening regime, presum- 
ably followed by a sharp peak stress corresponding to the 
unstable coalescence of cracks as the sample fails along 
a shear band. These results are consistent with what 
experimentalists observe. 

Using the exact differential equation that controls the 
temperature in the theory, it has been demonstrated that 
the temperature is becoming even more negative at lo- 
cahzation which means that the temperature is always 
finite at localization. Unfortunately, the exact value Tc 
of the temperature at localization is difficult to obtain 
because it is a result of integrating the differential equa- 
tion from the initial conditions. Although this could be 
done numerically, we have instead used an approximate 
value of Tc based on a non-interacting crack model. 

By far the most important signature of the localization 
transition is the divergence of the aspect ratio of the crack 
clusters. As reported, no definitive experimental work 
has yet been performed to test this prediction and we 
hope that experimentalists take this as a chahenge. 
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